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АВТОМАТИЗОВАНЕ УПРАВЛІННЯ ІМУННОЮ ВІДПОВІДДЮ НА 
ЗБУДНИКІВ ІНФЕКЦІЙНОГО ЗАХВОРЮВАННЯ 


Анотація. Проблеми, що виникли в результаті швидкого поширення нової коронавірусної інфекції 
СОМІР- 19 та неспроможності закладів охорони здоров'я своєчасно надати потрібну медичну допомогу великій 
кількості пацієнтів, засвідчили необхідність розробки нових підходів щодо діагностування та формування 
персоналізованих програм лікування. Важливим у цьому напрямку є розробка і впровадження спеціалізованих 
систем підтримки прийняття рішень на основі інтеграції експертних систем аналізу 1 інтерпретації медичних 
даних та інструментарію математичного моделювання. Для формування ефективних програм лікування 
необхідно забезпечити наявність надійного інструментарію прогнозування динаміки захворювання в умовах 
зовнішніх лікувальних впливів, різних механізмів захисту, просторових ефектів, змішаних інфекцій тощо. 
Метою даної роботи є узагальнення базової моделі вірусної інфекції для прогнозування комплексного 
керованого впливу на розвиток захворювання адсорбційної та імунної терапії з урахуванням дифузійного 
збурення та температурної реакції організму, що є важливим елементом спеціалізованих систем прийняття 
рішень щодо формування ефективних програм лікування. Для знаходження розв'язку вихідної модельної 
сингулярно збуреної задачі із запізненнями модернізовано процедуру поетапного чисельно-асимптотичного 
наближення розв'язків відповідної спеціальним чином побудованої послідовності задач без запізнення. На 
основі отриманих результатів комп'ютерного моделювання підтверджено, що введення адсорбентів у періоди 
інтенсивного зростання кількості вірусних елементів дає можливість достатньо ефективно зменшити швидкість 
їх розмноження, а, отже, забезпечує додатковий спосіб керування кількістю антигенів в організмі та 
недопущення надкритичного їх зростання. Підкреслено, що урахування особливостей дії адсорбційних 
препаратів є важливим для прийняття оптимізованих рішень щодо застосування адсорбційної терапії. 
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АСТОМАТЕФР СОКТКОЇ, ОК ТНЕ ІММОХЕ КЕ5ЗРОХ5Е ТО ТНЕ 
САО5АКТ58 ОЕ ІХЕКЕСТІОД5 РІЗ5КА5Е 


АБеігасі. Тре ргобіетя Ша агобе аз а ге5ції ої Ше гарі4 зргеад ої Фе пем/ согопауіги5 іпіесйоп СОМІР-19 
апа Фе іпарійсу ої пеаї Юр саге іп5цийопя (о ргоміде Ше песез5агу тедіса! саге іо а Іагеєе питег ої райепіз іпа 
йтеЇу таппег 5помед Ше пееад (о деуеіор пем арргоаспез іо дйїаєпозіяє апа Ше Гогтайоп ої ап іпаїмідца! (геайтепі 
ргоєгат. Пппрогіапі їп іє дїгесйоп 15 Ше деуеіортепі апа ітріетепіайоп ої зресіайгед десізіоп 5иррогі зузіетя 
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Їог ргедісійпо Ше дупатіс5 ої їБе д4ізеа5е іп Ше сопаїоп5 ої ехіегпаї пледіса! іпйЙиепсе8, уагіой5 ргоесйоп 
теспапізтя, зрацаї ейесіз, тіхей іпіесйопя екс. ТБе ригроз5е ої ія могК 15 (0 бепегайле Ше Базіс плоде! ої уїтаї 
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тодегпі7ед. Вазед оп Ше обіаїпеай гезиіз ої сотриїег подейпє, її маз сопіїгтеад Ша Ше іпігодисйоп ої ад5огрепіз 
Фитіпє регіодз8 ої іпіепзіуе гом ої Фе пипрег ої уїга! еІ|етепіз паКе5 іс ро55їбіе о 5иййісіепйу еНеспує!у гедисе 
Фе гаїе ої ШБеїг гергодисйоп, апа, ірегеїоге, ргоуідеє ап аддїйопа! мау ої сопігойтпя Ше пипег ої апійдепя іп Ше 
роду апа ргеуепйпе Шеїг зирегсгійса! єтгом'ФЮ. І 15 етарраз8ігед іа какіпє іїпіо ассоцпі Ше Теашгез ої ре асйоп ої 
ад5огрйїоп дгиєз 15 ітрогіапі ог паКіпє орійтігей десізіоп5 герагаїпе Ше цп5е ої адзогріїоп ШФегару. 


Кеуууогаз: іпіеспоця дїзеа5е тоаеі, ад5огріїоп Фегару, азутаріобс теїоайз. 


Вступ 

Виникнення та швидке поширення 
нової коронавірусної інфекції СОМІР-19, 
що було визнано надзвичайною ситуацією 
міжнародного значення та спричинило 
запровадження безпрецедентних заходів у 
всьому світі, 3 поміж іншого, стало 
свідченням того, наскільки високий рівень 


загрози з для людства шннесуть в собі 
інфекційні | захворювання. Поява у 
короткий термін великої кількості 


інфікованих СОУІД)-19 також вказала на 
неготовність наявних національних систем 
охорони здоров'я надати своєчасну та 
належного рівня медичну допомогу усім 
хворим. Цю проблему посилювало 1 те, що 
потужність реагування імунних систем 
різних пацієнтів на нову коронавірусну 
інфекцію виявилась досить різною, що 
зумовило появу 1 різних форм прояву цієї 
хвороби та можливість розвитку у них 
досить широкого спектру ускладнень. Такі 
особливості протікання СОМІ)-19 1 інших 
інфекційних 0 захворювань зумовлюють 
необхідність більш швидкого реагування 
на них, що включає, зокрема, як діагносту- 
вання вірусної інфекції на стадії (її 
розвитку в організмі, так і формування 
відповідної оптимізованої та 
персоналізованої | програми лікування. 
Одним із напрямків вирішення цієї 


проблеми є | розробка і широке 
впровадження у медичну практику 
спеціалізованих систем підтримки 


прийняття рішень на основі інтеграції 
експертних систем з аналізу і інтерпретації 
медичних даних та інструментарію 
математичного моделювання. 


Постановка проблеми 

Зазначимо, що описана в |Ї1| так 
звана | найпростіша модель вірусного 
захворювання забезпечує можливість для 
прогнозування загальних закономірностей 
протікання хвороби з урахуванням базових 
механізмів імунної відповіді гуморального 
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типу. Незважаючи на невелику кількість 
змінних, ця модель, з поміж іншого, дає 
можливість, залежно від ефективності 
імунної відповіді на вірусні елементи, 
виділити різні форми протікання хвороби: 
субклінічну, гостру, хронічну та летальну, 
що суттєво впливає на вибір стратегії 
лікування захворювання. А побудовані за 
загальною методологією моделі противі- 
русної та протибактеріальної імунної від- 
повіді (1| забезпечують ще й можливість 
прогнозувати розвиток відповідних видів 
захворювань з урахуванням особливостей 
як гуморального, так і клітинного типів 
імунного захисту організму. Поряд з цим, 
при формуванні ефективних програм ліку- 
вання з допомогою спеціалізованих систем 
підтримки прийняття рішень необхідно 
також забезпечити наявність інструмента- 
рію для прогнозування динаміки захво- 
рювання в умовах зовнішніх лікувальних 
впливів, інших механізмів захисту, зміша- 
них інфекцій, просторових ефектів тощо. 

Аналіз останніх досліджень |і 
публікацій 

Наведені в (|1| вже класичні моделі 
імунної відповіді створюють основу для 
розробки різного роду їх узагальнень і 
модифікацій |2|. Зокрема, в |3| наведена 
модифікація базової моделі інфекційного 
захворювання для побудови дискретного 
керування імунною відповіддю, а в |4| 
представлено підхід для побудови 
розподіленого керування, що забезпечує 
більш швидкий вихід на стаціонарний 
стан. В роботі |5| на основі спеціальним 
чином модифікованої моделі противірус- 
ної імунної відповіді досліджено різні 
сценарії імунної реакції організму на 
коронавірусну інфекцію  СОМІП2-І9 в 
умовах імунодефіциту та імунотерапії. 

В |6| представлено підхід для ураху- 
вання впливу на розвиток вірусної інфекції 
малого дифузійного розсіювання діючих 
чинників, а також показано, що урахуван- 
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ня цього ефекту призводить до очікувано- 
го зменшення модельної концентрації 
антигенів в епіцентрі зараження, а отже і 
до зниження прогнозної «гостроти» проті- 
кання хвороби. В |7) такий підхід розви- 
нуто для узагальнення моделі класичної 
противірусної імунної відповіді, а в |З, 9| 
класичну модель інфекційного захворюва- 
ння модифіковано для урахуванням дифу- 
зійних збурень відповідно в умовах фарма- 
котерапії та імунотерапії. Поряд з цим 


потрібно зазначити, що моделі 
інфекційних захворювань з дифузійними 
компонентами | на сьогодні є ще 
малопоширеними. Приклади таких 


моделей наведені в (10-14). 

Для прогнозування терапевтичних 
ефектів, що виникають внаслідок введення 
в організм лікарських препаратів (зокрема, 
за допомогою шприца чи крапельниці), в 
П5)| запропоновано модифікацію базової 
моделі шляхом введення додаткових ком- 


понентів, які дають можливість 
враховувати відповідні зосереджені 
впливи. А в П6) представлено 


узагальнення моделі противірусної імунної 
відповіді для комплексного урахування 
дифузійних збурень, температурної реакції 
організму, логістичної динаміки антигенів 
та зосереджених впливів. 

Поряд з цим зазначимо, що згідно 
засвідченої медичної практики застосуван- 
ня адсорбційних препаратів в умовах 
інфекційного захворювання дає можли- 
вість підвищити лікувальний ефект тради- 
ційних терапевтичних методів, а також 
забезпечує додатковий інструмент для 
контролю за розвитком хвороби. 


Мета дослідження 

Метою даної роботи є узагальнення 
моделі інфекційного захворювання з ура- 
хуванням дифузійних збурень, зосередже- 
них впливів та адсорбційної терапії для 
побудови ефективного управління імун- 
ною відповіддю на збудників вірусного 
захворювання. 


Виклад основного матеріалу 

Розвиток процесу інфекційного за- 
хворювання з урахуванням дифузійного 
збурення в умовах температурної реакції 


12 


організму та застосування адсорбційних 
препаратів для зручності викладу опишемо 
(для | простоти  викладок) в області 
СаО) -о«хчо, 1»0)| такою сингу- 
лярно збуреною системою нелінійних ди- 
ференціальних рівнянь із запізненням т: 


2 
ак, -(В(8)-укуУ єр ---- 1 
ОЇ дх 


ої ем 
2 

М тери 2 2 і 
ОЇ ді тд 


со -Бту ВУРИ-АУУ (т) 


з (1 


ж С 
-нс(С-С уре ро до 


ОК 
зр РЕСУР 
2 
Зерно. 
дх 
дт 2руп т 
-сУ- фер 
др -с ит Є зас? 
08 078 


т УЕ-ц (0-8) єр, 15 аа 


за умов 
Сх, сх), т(х, 0) ят'(х), 
0(х, 0) 07 (0, У Со) «У"(х, 1), 


Е(хї)з Е"(х,1), -т«ї «0, (2) 
У (ж, г, 0)- МУ" (хг), НОЙ, 


У (х, Кут) ек У (ху), 
де У(х,ї), К(х,і), С(х,ї), т(х,ї), 0(х,1) - 
відповідно концентрації антигенів, специ- 
фічних їм імунних агентів, імунологічних 
клітин, значення відносної характеристики 
ураження органу-мішені та температури у 
точці х в момент часу І; М/(х,г,І) - кон- 
центрації антигенів в мікропорах частинок 
адсорбенту (для простоти викладок вважа- 
тимемо, що вони мають форму куль радіу- 
са К з центром в точці х); к»0 -- 
константа адсорбційної рівноваги; В(0)- 
«в/ч, (0-6) - залежний від темпе- 
ратури організму темп розмноження анти- 
генів, В, -сопіт»0; у - коефіцієнт, що ви- 
значає результат взаємодії антигенів З 
їмунними агентами; г - запізнення в часі; 
ше - величина, обернена тривалості життя 
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плазматичних клітин; 


а(8) з а (Ічна, (0-8"))  - залежний від 
температури організму коефіцієнт 
стимулювання імунної системи, 


а, -сопії»0; С - концентрація імуноло- 


гічних клітин у здоровому організмі; р - 
швидкість виробництва власних імунних 
агентів однією імунологічною клітиною; г; 
- величина, обернена тривалості існування 
імунних агентів; 7) - витрати імунних 
агентів на нейтралізацію одного антигену; 
с -- темп ураження клітин органу-мішені; 
ИШп - швидкість відновлення органу-міше- 

; (9) - розподіл температури в здоро- 
вому організмі; 


 Г0, УКА«(УБ), 
"Дао, УК2(УЕ); 


Т? 


(УК) - порогове значення УК - комплек- 
сів, при якому підвищення температури ще 
не стимулюється (тут прийнята гіпотеза, 
що температура тіла залежить від концент- 
рації УК - комплексів), ау зсопіт»0; єр, 
єр", єр," єр, - коефіцієнти дифу- 
зійного розсіювання відповідних анти- 
генів, імунних агентів, імунологічних та 
уражених клітин, теплопровідності, є р" 
- коефіцієнт дифузійного розсіювання 
антигенів в мікропорах адсорбенту, єр" - 
коефіцієнт, що характеризує вплив 
дифузійного розсіювання антигенів у 
мікропорах адсорбенту на їх дифузійне 
розсіювання в середовищі організму, Є 

параметр, що характеризує малий вплив 
дифузійних та адсорбційних компонентів 
на розвиток процесу у порівнянні з іншими 


його складовими; С"(х), т'(х), 67(х), 
у?(х,ї), Е"(ої), У/"(хуг) - обмежені 
достатньо гладкі функції. Функція 5(т) 


забезпечує урахування ефекту зниження 
інтенсивності розмноження імунологічних 


клітин при значному ураженні органу- 
мішені, функції с (хг), и" (01) 
призначені для урахування, зокрема, 


зосереджених змін концентрацій антигенів 
та імунних агентів (15, 16). Функцію 
У 

и" (0,1), забезпечує 


яка «запуск» 


механізму поглинання антигенів частинка- 
ми адсорбенту у визначені моменти часу, 
представимо у внрияль 


и" (те о, 


Для знаходження розв'язку задачі 
(1)-(2) із запізненням, як і в |6-9, 15, 16, 
скористаємось відповідною  покроковою 
процедурою, замінивши її послідовністю 
задач без запізнення на кожному з проміж- 
ків 5ссІХ(5Ч1)с, 50,1... із забезпечен- 
ням необхідної гладкості цього розв'язку 
при іт, 27,..! 


ет (и, 


ре (и) 


2 


бу б у 0 Ха 
І чо (В0 УК У те 


се 
дт рек 


Р го (Зо ато 


др об 


З) 
2: 
-С'узе рі Су 


"Ше 


2 


дт д'т 


98. запро 


б "за то Нт (8 


С уСьзт)е 


туру 557), 0 (57) 26, (357); 


Су 557), т, (х, зт) 2 


сю 51) - У у 5 87), йон (4) 
ЕС 5), о (г 8т) ЗМ СХ ЗТ), 
а050,7) 
УУзбьВгеку,, (г), стано 0, 
де С (х,б)єС"(9, | т (хбувт'(ю), 


0. (2, 0)-0"(х), У (0) -У"(х,0), 

Е (о0)-к"(х0), || У (хугОу- М" (хг, 

ню по о (з 
РОМЕО г з ЛЕО 2 сі А і РРО 


Зважаючи на те, що дифузійне розсі- 
ювання діючих чинників процесу є малим, 
для наближення розв'язків часткових син- 
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гулярно збурених задач (3)-(4) на кожному 
з вказаних вище часових проміжків засто- 
суємо  чисельно-асимптотичний метод, 
формально представивши їх у вигляді 
асимптотичних рядів |6-9, 15, 16. 
Автоматизоване формування ефек- 
тивної програми лікування інфекційного 
захворювання передбачає, зокрема, роз- 
роблення такого управління імунною від- 
повіддю організму, яке забезпечуватиме 
рівень гостроти протікання хвороби та 
темп відновлення ураженого  органу- 
мішені у наперед визначених допустимих 
межах. Відшукання такого управління 
зумовлює наявність інструментарію надій- 
ного прогнозування розвитку хвороби за 
різних ситуаційних станів та умов засто- 
сування різного роду керуючих зовнішніх 
впливів. Наявність в організмі різних 
захисних механізмів ускладнює класичну 
постановку задачі управління відповідною 
комплексною реакцією на виявлені вірусні 
елементи з однозначно обраним критерієм 
оптимальності. Крім того, в реальних ситу- 
аціях досить складно з потрібною точністю 
спрогнозувати силу імунної відповіді в 
умовах застосування різних видів ліку- 
вальних препаратів, а також різних спосо- 
бів їх введення. Таким чином, для корегу- 
вання терапевтичного впливу на розвиток 
інфекційної хвороби виникає потреба 
періодично в процесі лікування оцінювати 
рівень ураження організму та стан імунної 
системи. В цих умовах контроль та корек- 
цію динаміки інфекційного захворювання 
раціонально здійснювати із застосуванням 
дискретного адаптивного управління імун- 
ною відповіддю |17). Згідно цього підходу 
спочатку визначається траєкторія  «ба- 
жаної» імунної відповіді у вигляді множи- 
ни точок У(г)-У,, і-1,М. Далі з набору 
кусково-сталих на відрізку І0,7| функцій 
0 АЧи(): и)зи,, ЄЇ0, В), тег, ,,1,),1-1,М, 
и(Т)зи, |) За результатами інтегрування 
вихідної модельної задачі з управління на 
проміжках |/. ,,1,|, і-1,.,М при фіксова- 


них значеннях и,, визначають управління 
и(т), яке забезпечує мінімізацію функції 
Ф(и, |)- Ми, 1)- У, . 

У результаті застосування адсорбцій- 
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них препаратів в моделі (1)-(2) «запуск» 
механізму поглинання антигенів 
адсорбентами описується функцією и" для 
управління динамікою вірусного 
захворювання, при якій очікується 
зниження кількості збудників цієї інфекції. 
А, отже, для знешкодження меншої 
кількості антигенів необхідна буде і менша 
кількість | антитіл, | які | виробляються 
імунною системою чи вводяться шляхом 
ін'єкцій (введення в організм донорських 
антитіл у моделі (1)-(2) описується ще 
однією функцію управління и" ). Таким 
чином, відповідне дискретне адаптивне 
управління імунною відповіддю на збудни- 
ків вірусної інфекції в умовах адсорбційної 
терапії та температурної реакції організму 
представимо у вигляді такого алгоритму. 
На початку почергово для кожного про- 
міжку |? ,,1Ї, 1-1,..,М знаходимо управлі- 


У Б | здоб 
ння и, яке, за умови и ,-0, мінімізує 

: 372 г У Е з 
функцію  Ф(и: ,и,,|)-У(ицц)-У . У 


випадку, коли за такого управління на 
цьому проміжку досягається рівність 


177 ГЕ 177, ж 
Ф(и, и, |) - У (ці ,0,)-У «0, переходимо 
до наступного проміжку. В іншому ж 
випадку - фіксуємо знайдене значення 

У/ А й 
из; та переходимо до відшукання зна- 
Е зло є Р 

чення и; ,, яке забезпечує мінімізацію фун- 
У Б л ГЕ б 
Фан, )2 У (Цой) - 
цьому, природно вважати, що функції и" 
БЕ 177, У 
та и є Охи «их» 


Б В 
0 А ей з 


кції При 


обмеженими: 


Числові експерименти 

Здатність адсорбційних речовин по- 
глинати вірусні елементи забезпечує мож- 
ливість здійснювати додатковий коригую- 
чий вплив на розвиток вірусної інфекції, 
що є досить важливим, зокрема, у періоди 
загострення хвороби. Проведення число- 
вих експериментів були зосереджені на до- 
слідженні прогнозного впливу адсорбції на 
модельну динаміку інфекційного захворю- 
вання. 

На рис. І представлено прогнозні ди- 
наміки кількості антигенів в епіцентрі за- 
раження у випадках, коли: ефект дифузій- 
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Рис. 1. Динаміка антигенів в епіцентрі зараження при Є)" «0.00 і ЄЮ" «0.00 (крива 1); ЄДІ" «0.05 і 
єр «0.00 (крива 2); ЄДІ" «0.05 і є Р" «0.05 (крива 3) 


ного розсіювання діючих чинників не вра- 
ховується та відсутня адсорбція (крива 1); 
дифузійне розсіювання діючих чинників 
враховується, але адсорбція відсутня (кри- 
ва 2); ефект дифузійного розсіювання вра- 
ховується та має місце адсорбційне погли- 
нання антигенів, яке починається у момен- 
ти часу І, 79» Їу, -8 і їу, 7-12 (крива 3). Як 


і очікувалось, в останньому модельному 
випадку спостерігається найнижче макси- 
мальне значення прогнозної кількості ан- 
тигенів, внаслідок чого забезпечено | від- 
повідне зниження критичного рівня загост- 
рення інфекційного захворювання. Резуль- 
тати числових експериментів також засвід- 
чили те, що ефективність застосування ад- 
сорбентів для стримування росту популяції 
антигенів залежить від моменту часу їх 
введення. Більш ефективний вплив на 
швидкість розмноження антигенів 
отримано у випадках, коли введення 
адсорбентів з здійснювалось у періоди 
збільшення їх кількості. Урахування такої 
особливості дії адсорбційних препаратів є 
важливим для прийняття оптимізованих 
рішень щодо застосування адсорбційної 
терапії. 


Висновки 

У роботі узагальнено модель вірусної 
інфекції для прогнозування керованого 
впливу адсорбційної дії та імунотерапії на 
розвиток захворювання в умовах дифу- 
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зійних збурень та температурної реакції 
організму, що є важливим елементом спе- 
ціалізованих систем прийняття рішень. 
Для знаходження розв'язку вихідної мо- 
дельної сингулярно збуреної задачі із за- 
пізненнями модернізовано поетапну про- 
цедуру чисельно-асимптотичного набли- 
ження розв'язків відповідної послідовності 
задач без запізнення. 

Результати числових експериментів 
демонструють додаткове зниження прог- 
нозної кількості антигенів в результаті зас- 
тосування адсорбційної терапії. Підтверд- 
жено, що введення адсорбентів у періоди 
інтенсивного зростання кількості вірусних 
елементів забезпечує достатньо ефективне 
зменшення швидкості їх розмноження, а, 
отже, і відповідне зниження надкритичної 
кількості антигенів та більш легке проті- 
кання хвороби. Наявність інструментарію 
для прогнозування динаміки інфекційного 
захворювання з урахуванням керованих 
впливів адсорбційної та імунної терапії в 
умовах температурної реакції організму та 
дифузійних збурень у спеціалізованих сис- 
темах прийняття рішень дає можливість 
формувати відповідним чином оптимізо- 
вані програми лікування. 

У перспективі є розвиток представле- 
ного підходу для урахування різних меха- 
нізмів дії інших типів сорбційних препа- 
ратів, конвекції, умов змішаних інфекцій 
та керованих зосереджених впливів фарма- 
котерапії 1 біостимуляції. 
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